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Abstract: Chandra [1] made an attempt to show that the work of Dwivedi and 

Srivastava [2] (hereinafter DS) can be investigated even analytically and their results 

are erroneous. Dwivedi and Srivastava [3] picked up some values of Chandra [1] and 

tried to show that they are not physically acceptable. Some results of Chandra 

[1] are not physically acceptable, as these are the outcome of the wrong approach 

of DS. However, the results are numerically correct whereas the results of DS are 

i-C ■ numerically wrong. 
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c3 ! 1 Introduction 



In order to show that the work of Dwivedi and Srivastava [2] (hereinafter DS) can be 
investigated even analytically, Chandra [1] used the condition ujv >> v\ for the case 
of rj = and ur) « v\ for the case of v = 0. His also proved that the results of DS 



in 



o 
o 



are wrong. Dwivedi and Srivastava [3] should not misguide the scientific community 
by wrongly interpreting these conditions, as the numerical calculations can also 
prove that the results of DS are wrong. Moreover, as a consequence of the wrong 
approach of DS, some numerical results obtained by Chandra [1] are not physically 
acceptable. In the present work we have shown numerically that the results of DS 
are wrong. 
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2 Formulation 



The dispersion relation derived by Pekiinlu et al. [4] for the MHD equations can be 
expressed as [1, 2, 4] 

z^A; 4 + [v\ - iu(y + r])]k 2 - u? = (1) 

Consideration of an atmosphere with (i) v = (no viscosity) or (ii) with r\ = (no 
magnetic diffusivity) may not be a physically acceptable situation. However, DS 
discussed these two cases numerically, where they got wrong results. Let us now 
investigate these two cases numerically. For either of the cases, the first term vr\k 
is zero and we have 

v\ - + r])]k 2 - u? = (2) 



It is matter of surprise that Dwivedi and Srivastava [3] are asking for solving equation 
(refeqO) numerically. For real values of lo, the wave number k can be assumed to 



have complex values. For k = k r + iki, we can derive the values of k r and ki for the 
said two cases as the following: 

(i) There is viscosity only . That is, there is no magnetic diffusivity (77 = 0). For 
this case, equation (2) gives [1, 2] 



kf-P k 2 r -Q 2 = 



and 



— 



(3) 



where 



P 
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Form equation (3), we get 



k 



2„,2 



LU V 



and 



Q 



2{v\ + u 2 u 2 ) 



r 2(v\ + u 2 u 2 ) 



1 + 



yJl + (u}v/v 2 A )< 



(4) 



(ii) There is magnetic diffusivity only. That is, there is no viscosity [y = 0). For 
this case, equation (2) gives [1, 2] 



fc r 4 - P'k 2 r -Q' 2 = 



and 
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(5) 



where 



P' 



v\ + UJ 2 f} 2 



Form equation (5), we get 



2 2 



and 



2(«J+«V) 



2 \2 



(6) 



3 Results 



In order to get numerical values for k r and ki, for the two cases, we require values 
of physical parameters. The values of physical parameters are provided by Chandra 
[1] and some of them are given in Table 1. In order to find out the source of error 
in the work of DS, we asked them several times to provide us the values of physical 
parameters used in their work, but they never showed the courtesy to reply. For 
r = 10~ 4 s, the calculations are made and the results are given in Table 1. The 
damping length scale D and the wavelength A of the wave are 



D 



2tt 

ki 



and 



A = 



2tt 



The results for the case r\ = show that the damping length scale and the wavelength 
are of the same magnitude, which is of the order of 10 4 m. For the case v = 0, the 
damping length scale comes out of the order of 10 11 m.This physically unacceptable 
result is obviously a consequence of the wrong approach of DS. Comparison of the 
present results with those of DS show that their results are wrong. Even the nature 
of the variation of D is quite different. 



Table 1: Values of k r and ki for the two cases. 
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"Distance from the center of sun in i? Q 
fc Coemcicnt of viscosity in 10 11 m 2 s _1 ; 
c Magnetic diffusivity in m 2 s _1 
d Alfven velocity in 10 6 m s _1 
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